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1. Introduction 

I describe a strictly local fermion Dirac operator &{A) with an exact chiral symmetry mani- 
fested in anti-commutation with 75, i.e. 75^ = —^75. By strictly local I mean that only nearest 
neighbor terms appear in the fermion action. The operator describes two fermion flavors, the min- 
imum required for chiral symmetry to exist. I develop this action as a linear combination of two 
"naive" fermion actions, following a line of reasoning similar to that presented by Borici [[I]]. 

The theory is not symmetric under the full hyper-cubic group, but the subgroup thereof that 
preserves one fixed direction up to a sign. These symmetries include transformations of both even 
and odd parity. On renormalization, interactions can introduce a lattice anisotropy at finite cutoff. 



2. Doubling and chiral symmetry 

Spontaneously broken chiral symmetry is fundamental to our understanding of low energy 
hadronic physics. Pions are elegantly described as quantum mechanical waves propagating through 
a background quark condensate. In addition, chiral symmetry provides powerful tools enabling 
extrapolations to the physical quark masses from the heavier values currently practical in lattice 
gauge simulations. 

These issues are deeply entangled with quantum mechanical anomalies that eliminate one 
symmetry of the classical Lagrangian. With Nf flavors of quark, the naive axial U {Nf) is reduced 
to SU (Nf). With only one flavor, all chiral symmetry is removed; so, multiple flavors are necessary 
for chiral symmetry to be relevant. Nielsen and Ninomiya |§] have given a formal topological 
argument that any lattice action with chiral symmetry must describe at least two flavors. 

If one ignores the anomaly and writes a simple lattice action that is chirally symmetric, some- 
thing must go wrong. The usual result is that the fermion field describes multiple species, and 
the extra species cancel the anomalies. The most naive discretization, which will play an essential 
role below, involves 16 species in four dimensions. Staggered fermions divide out an exact SU (4) 
symmetry of the naive formulation to reduce the multiplicity to four [||]. The Wilson fermion 
approach successfully removes all doublers at the expense of breaking all chiral symmetries. El- 
egant newer approaches based on perfect actions [Q], domain walls in five dimensions [^], or the 
overlap operator do maintain much of the desired chiral symmetry with arbitrary Nf, but involve 
computationally expensive non-local actions. Also, these approaches tend to obscure the anomaly; 
for example, with the overlap one introduces a gauge field dependent matrix 75, the trace of which 
gives the winding number of a given gauge configuration. 

The Nielsen-Ninomiya theorem requires any chirally symmetric lattice action to describe at 
least two species. Actions which satisfy the minimal doubling of just two have been known for 
some time ||8|, ^], and have recently stimulated new interest [|], |lO, II, 12]. There are a variety of 



compelling reasons to revisit these actions. First is the failure of the rooting procedure popularly 
used to reduce the doubling of staggered fermions [|T^]. Second is the lack of an exact chiral 
symmetry for Wilson fermions, complicating extrapolations to physical fermion masses. And third 
is the severe computational demands of the domain-wall and overlap approaches. 
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Figure 1: Representing momentum space as a product of toroids, a massless Dirac equation can be obtained 
from naive fermions by expanding about any point where the components of momentum are either or ;r. 
This gives rise to the 16 doublers of the approach. 

Here I construct a minimally doubled fermion action closely following Borici's formalism 
I use a linear combination of two unitarily equivalent naive fermion actions. This combination will 
be crafted so that only two of the original 16 doublers for each action survive. 

3. Naive fermions 

The so called "naive fermion" approach plays a crucial role in the following; so, I begin by 
reviewing this action. Work with a conventional hyper-cubic lattice with gauge fields implemented 
as group elements on the links. The details of the gauge part of the action play no role here. When a 
fermion hops in a forward direction jx between neighboring sites, it picks up a factor of YfiU . Here 
7^ is the usual Hermitean Dirac matrix while U represents the gauge field on the corresponding 
link. For the reverse, or backward, hop on the same link, the contribution is —y^U'^. The minus 
sign makes the fermion operator an anti-Hermitean matrix when viewed as a matrix in the direct 
product of the site space with the internal symmetry and spinor spaces. I work here with the 
convention of a unit hopping parameter; for the massless case any other hopping parameter can be 
scaled away with a redefinition of the fields. I also work in lattice units so that aU site coordinates 
are integers and the lattice spacing does not appear explicitly. 

As is well known, this action describes the physics of 16 fermion species, frequently referred 
to as doublers. The Dirac operator, D, anti-commutes with 75. This represents an exact chiral 
symmetry. Because the doublers use different effective gamma matrices, half of them rotate in 
each direction under a chiral rotation. Thus this is a non-singlet chiral symmetry. 

In the free field limit where all the hnk matrices are the identity, this theory is easily solved 
in momentum space. The Dirac operator factors into independent pieces for each momentum p, 
taking the form 

D(/7)=2/£7^sin(/7^). (3.1) 

Expanding about small momentum gives the usual Dirac behavior D{p) ~ 2/^^ y^iP^i- The dou- 
blers appear on expanding not about zero momentum but around points where some components of 
p are approximately n. Thus, in addition to the excitations around p^ = 0, there are 15 other points 
in momentum space where the action is small. Visualizing momentum space as a direct product of 
two toroids, the zeros occur as sketched in Fig. ([T]). 
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Figure 2: The points in momentum space furthest from those giving the Dirac equation are located at 
Pfi = ±7r/2. Our second naive fermion action will have its zeros at these points. 

4. A unitary transformation 

I now do a transformation on the above action to generate a superficially different but physi- 
cally equivalent naive fermion action. Begin by considering maximally distant momenta from the 
zeros of D. There are 16 such points, occurring whenever the components of the momentum satisfy 
= -iin/l. These points are sketched in Fig. for the x,y sub-torus. Arbitrarily select one of 
these points; here I consider = +n/2 for all pL. Here the original action becomes 



where I define the quantity 



Z)(p^ = 7r/2)=2/£7^=4/r 



1 . 

r= -(71 + 72 + 75 + 74)- 



This is a unitary, Hermitean, and traceless four by four matrix in spinor space. 
Now consider a unitary transformation on the original fields 

^ff = g-in{xi+X2+XT,+Xii,)/2 Y Y 

-j^ _ gi7C{xi+X2+X3+X4)/2 -j^ Y_ 



(4.1) 



(4.2) 



(4.3) 
(4.4) 



Here the are the integer coordinates of the lattice. The phases move the zeros in momentum 
space from p^ =0,71 to the maximally distant points p^ = ±7r/2. The factors of F modify the 
gamma matrices for the new action to 

/^=F7^F. (4.5) 
Note that I can construct F either from the original 7^ or the new 

r = ^(71 + 72 + 73 + 74) = ^(71 + /i + T^s + T^) = r'. (4.6) 
In the free field limit, momentum space again diagonalizes the new action, call it D, 

D{p)=2i'£y^sm{n/2-p^). (4.7) 

Of course since all I have done is a unitary transformation, the Dirac operators D and D are physi- 
cally equivalent. They also satisfy a complementarity manifested in 



D{p^ = K/2)=D{p^=0)=4ir. 
This property is crucial to the following construction. 



(4.8) 
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5. The minimally doubled action 

I now construct the final action from a linear combination of these two equivalent theories 

& = D + D-AiT. (5.1) 
To see how this works, consider the free theory limit in momentum space 

i^(p) =2/£(7^sin(p^) + /^sin(7r/2-p^)) -4/r. (5.2) 

Going to the Fermi point of the original theory at ~ 0, then the 4/r term cancels D. On the 
other hand, at the Fermi point of D occurring at ~ ;r/2, the 4/r term cancels D. The remarkable 
feature of this combination is that only these two zeros of ^(p) remain. I give a more detailed 
proof in Appendix A, but basically at the other zeros of D, D — 4/r is large and at the other zeros 
of D, D — 4/r is large. Fig. (^) shows the two remaining Fermi points. 

Note that each term in Eq. ( p7T[ ) anti-commutes with 75, maintaining the exact chiral symmetry 
of the naive action. Thus a finite chiral rotation gives 

e'Oy^^e'^rs (5.3) 
The construction uses different gamma matrices for the two species. A particular consequence is 



= TysT = —75. Under the chiral rotation of Eq. (5.3), the two species behave oppositely. The 
symmetry is of a non-singlet nature, as expected. 



6. Space time symmetries 

The above construction utilizes the matrix r = jY.^Y^. Had I selected another of the max- 
imally distant points from the zeros of the initial fermion operator, then this relation would be 
modified with minus signs for some directions. This would give an equivalent theory, but any such 
choice still involves picking a particular diagonal axis of the fundamental hyper-cubic lattice as 
special. Here I have chosen the positive major diagonal. For comparison, in Refs. [jS], ^] the special 
direction is not a diagonal but rather X4, as discussed further in Appendix B. 

Having picked a specific diagonal of the fundamental hypercubes, the action is not symmetric 
under the full hyper-cubic group, but only the subgroup that leaves this direction invariant. This 



5 



Local chiral fermions 



Michael Creutz 



subgroup includes Z3 rotations that cycle any three positive axes. For example, I might want to 
cyclicly permute the xi,X2,xj, axes. To get the gamma matrices to transform appropriately under 
such rotations, introduce the matrix (in spinor space) 

y = exp ^^^(ai2 + a23 + a3i)^ . (6.1) 

Here I define [7^,7v]+ = 2ia^v The combination V^y^V cyclicly permutes the first three gamma 
matrices. Since [V,r] = 0, this matrix, along with a corresponding rotation of the gauge fields, 
realizes a symmetry of the action. Note that = —1, indicating that a rotation by an angle of 2 TT 
gives a minus sign, as expected for fermions. Combining such rotations using other axes generates 
the 12 element tetrahedral subgroup of the hyper-cubic group. 

The above rotations generate positive parity permutations of the axes. Introducing negative 
parity permutations increases the subgroup to 24 elements. For example, the matrix 

V = ^ ( 1 + /ais) ( 1 + /a2i ) ( 1 + /C752) (6.2) 

generates the fermionic part of a rotation that exchanges the xi and X2 axes. This transformation 
flips the sign of 75, i.e. V^/sV = —75, emphasizing that it is a negative parity transformation. 

In this formalism the natural direction to represent time is the main diagonal ei + 62 + e^, + 64. 
The combination of time reversal and parity can then be chosen to flip the sign of all axes. A 
unitarity transformation similar to that relating D and D restores the action to its original form, 
increasing the symmetry group to 48 elements. 

Note that charge conjugation symmetry is trivial in this formulation, just being particle hole 
symmetry. Both the operators & and the Hermitean combination 75 have their eigenvalues in 
opposite sign pairs. 

Because of the special treatment of the main diagonal, the effects of interactions can distort 
lengths along this direction. At each of the Fermi points, this distortion is associated with the di- 
mension five continuum operator nj/FV^Y- Interactions at finite lattice spacing can result in the 
gluons and fermions not having the same speed of light. This can be corrected with a renormal- 
ization of the ^rrn in the action, as emphasized in Ref. [pT|]. Nevertheless, the zeros are 
stable under this distortion because they involve a topologically non-trivial mapping of surfaces of 
constant action ^ [lO|]. 

Appendix A: Proof that there are only two zeros of ^(p) 

From the definition of F it is elementary to show that = F — 7^ and the properties TrF7^ = 
TrF]^ = 2. Using trivial trigonometric identities one can deduce 

Tr (7^ - rv)^{p) ~ sin(;;^ - 7i/4) - sm{py - n/4). (6.3) 

This implies that at any zero cos(/7^ — 7i/4) = ±cos{py — 7i/4). But a zero also requires 

Tr r&ip) = ^ £cos(7?^ - 71/4) = 2^2 > 2. (6.4) 

Thus all these cosines must be positive. Therefore all components of are equal and either or 
7t/2, the two Fermi points of interest. 



6 



Local chiml fermions 



Michael Creutz 



Appendix B: Comparison with actions from Karsten and Wilczek 

Two other minimally doubled chiral fermion actions have been presented in Refs. [jSj, These 
actions are in fact equivalent to each other under a unitarity transformation y/ e^^^^'^l^\\r. For the 
free case, that action can be written 

4 3 

^= £ 7AiSinb^j) + 74£(l-cos(/7i))- (6-5) 

The last term removes all doublers from the naive action except those at = and = 0, n. 

The main difference from the action presented here is that X4 is now chosen as the special 
direction. The on site term is proportional to 74 instead of Y. As before the second species with 
momentum around P4 = 7i uses different effective gamma matrices, f = Y and )^ = — 74. As with 
the action from the main text, this gives ]^ = —75. Again the chiral symmetry is flavored. The 
symmetry of this system is the subgroup of the hyper-cubic group that preserves the fourth axis up 
to a sign. 
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